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^ : ABSTRACT 
< 

0\ ' The recent Chandra image of the Crab nebula shows a striking, axisymmetric 

^ ' polar jet. It is shown that jets are formed in axisymmetric, magnetized pulsar 
winds and that the jet luminosity scales relative to the total as {'joO'eg)~'^^^ , where 
(Teg is the ratio of Poynting flux to particle kinetic energy output at the equator 

^ ■ at the base of the flow and 70 the initial Lorentz factor of the flow. The results 

Tj- ■ are applied to the image of the Crab nebula, and the limit is set for the Crab 

^ ' pulsar of (Teg < 10^. It is argued that conventional pulsar theory needs to be 

O ■ reexamined in light of these limits. 

^ I Subject headings: ISM:jets and outflows — MHD — pulsars:general — stars:winds, 

?-j I outflows — supernova remnants 



1. Introduction 

Magnetic collimation of astrophysical jets was proposed over 20 years ago (Blandford 
1976, Lovelace 1976, Benford 1978). Since then, however, there has been considerable debate 
as to how well it works. Heyvaerts and Norman (1989) and Chiueh, Li, and Begelman (1991) 
proved that given enough time, all rotating, axisymmetric magnetized flows are coUimated 
by the toroidal magnetic field that is inscribed by the rotation. On the other hand, it 
was argued (Eichler 1993, Begelman & Li 1994, Tomimatsu 1994, Bogovalov 1997, 2001, 
Beskin, Kuznetsova & Rafikov 1998, Chiueh, Li & Begelman 1998) that the collimation is so 
slow (logarithmic) that in practice it could take an exponentially large radius to coUimate 
appreciably, so that all but a small core near the axis would remain uncollimated. Moreover, 
if the magnetic field dominates the flow, collimation would come at the expense of kink 
instability. The slowness of the collimation and/or its questionable stability was felt by 
some to make it unattractive. 
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For relativistic outflows with Lorentz factor 7^1, the problem is particularly severe. 
Firstly, the radius of curvature of the flow lines scales as 1/7^, so that for large enough 7, 
the radius of curvature is not only exponentially large, but exponential large with a huge 
exponent. Secondly, even if the flow started at modest 7, if it were magnetized, even a small 
degree of collimation would imply acceleration (Begelman & Li 1994) to larger 7, which 
places even more severe limits on collimation. 

However, the remarkable image of the Crab Nebula in the X-rays recently obtained by 
the CHANDRA Observatory (Weisskopf et al. 2000) convincingly shows a very well defined 
polar jet. It appears to be a) very well collimated but b) only a small part (~ 10^^ — 10^^ by 
visual inspection) of the total energy, which is mostly expelled in an apparently equatorial 
disk. This vindicates both the contentions that a) magnetic collimation can work and that 
b) it only works on the inner core of an ultra-relativistic outflow. Stability worries may 
be minimized if most of the jet has been collimated only within its most recent dynamical 
timescale. 

The question is now whether theory can accurately account for the proportions of the 
polar jet relative to the total pulsar wind, whether the predicted Lorentz factor of the jet is 
consistent with observational limits on relativistic beaming, and whether the implications of 
the pulsar wind observations are consistent with the theory of pulsar magnetospheres. 

In this paper we consider these questions. We show that the conventional view of the 
Crab pulsar would allow too miniscule a jet to be consistent with the observed one. We 
then show that relaxing the usual beliefs about the Crab pulsar and assuming only that 
the pair luminosity must be less than the gamma ray luminosity allows a jet fraction at 

most 10^^7q just barely consistent with the observations even for 70 ~ 1. Here 70 is 
the Lorentz factor at the base of the flow. In the following sections, we show that the jet 
luminosity scales relative to the total as (7oO"e(j)^^^^! where a^q is the ratio of Poynting flux 
to particle kinetic energy output at the equator at the base of the flow. We note in the 
discussion section the implications of this for theories of the Crab pulsar. 

2. Basic equations 

We flrst review the basic theory of relativistic, magnetized winds. At large distances 
(in fact beyond the fast magnetosonic point) one can consider the flow as a purely poloidal 
{v^ = O) flow in a purely toroidal magnetic fleld (cf. Li 1996). Indeed conservation of the 
angular momentum implies that the azimuthal velocity goes to zero with radius while the 
contribution of the poloidal fleld to the force balance becomes eventually negligible because 
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the poloidal field falls off faster than the toroidal one. In ultrarelativistic flows, hoop stress 
and electrical force nearly cancel each other so that the poloidal field pressure, while much 
less than either term in magnitude, is negligible only far enough from the light cylinder to 
be smaller than their difference. However, this distance still lies well within the termination 
shock, and the logarithm of the ratio between the two will be shown to play a role. The 
momentum conservation equation is given by 

p(7V- V)7V = peE+j X B, (1) 

where p is the proper mass density, v the plasma velocity, E and B are the electric and 
magnetic field strengths, correspondingly and Pe and j are the charge and current densities, 
correspondingly. Under the above assumptions one can write 

V = -yl; B = Be^- E = -v x B = -vBt- (2) 

where 1 is the longitudinal unit vector along the magnetic surface, is the azimuthal unit 
vector and t = 1 x e,^ is the unit vector in the transverse direction. Ee take the speed of hght 
to be unity. With the Maxwell equations, we can present the current and charge densities 
in the form ^ 

47rj = V X B = -V(r5) xe^; (3) 

47rpe = V • [B X v] = -^t ■ V{rB) -Bt-Wv- vB{e^-[V x 1]), 

where r is the cylindrical radius. 

The cross-field equation may be obtained by taking the dot product of Eq.(l) with t 
and using the above expressions for pe and j. After simple algebra, we obtain 



p(f - l)t ■ (1 ■ V)l = \v'B\e^.[V X 1]) - -t • V— I , 



(4) 



where 7 = (1 — f ^) is the flow Lorentz factor. Taking into account the definition of the 
curvature radius of the flux line 

1 



— = t • (1 • V)l = -t • [1 X [V X 1]] = -e^-[V X 1], 

^ + - =-t-V(ln— ), (5) 



Rc 

one can write finally 



Rc \ a J 7 
where the ratio of the Poynting flux to the kinetic energy flux is 



_ EB _ B^ 
Anpv^^ 47rp7^ 
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One can see from Eq.(5) that in an ultrarelativistic wind, 7^1, the flow hnes are nearly 
straight and coUimation is possible, if possible at all, only near the axis (Tomimatsu 1994, 
Chiueh et al. 1998, Bogovalov 2001). 

The Lorentz factor of the flow should be determined from the energy conservation along 
the flux tube. This equation may be obtained by taking the dot product of the momentum 
equation (1) with the longitudinal vector 1. However one can instead simply write down the 
energy flux per unit poloidal flux 

Pl^v + ^ 

— ^ = const. 

Bp 

Taking into account conservation of the mass flow within the field line tube, we define a 
quantity 77 (^) that is conserved along each surface ^ by 

p^;7 = r]{^)Bp. (7) 

Making use of this equation and the relation 

E = rQSp, (8) 

one can write the energy equation as 

Here the poloidal flux function is defined from the poloidal field strength by 

Bp = - V* X e^. 
r 

It follows from Eq.(9) that plasma accelerates if rlV^'l decreases (Begelman & Li 1994). In 
both radial flows and cylindrical ones the Lorentz factor remains constant. (Note that in 
the case of interest one need not introduce the poloidal field at all but instead define ^ as a 
stream function (cf. Contopoulos 1995). We retained the standard definitions to provide a 
simple comparison with the previous works.) 

The asymptotic transfield and energy equations were found by expanding the full set of 
MHD equations in 1/r by Begelman & Li (1994), Tomimatsu (1994), Bogovalov (1997, 2001), 
Chiueh et al (1998), Bogovalov & Tsinganos (1999). One can show that their equations are 
equivalent to Eqs.(5, 9). Examining the neglected terms one can see that these equations are 
valid at r 3> 7-Rl where Rl = c/Vtis the light cylinder radius and the angular velocity of 
the pulsar. This condition arises because the neglected poloidal field stress (oc oc [r/RiY) 
should be compared in the relativistic case not with the hoop stress (oc B'^ oc (r/RLY) but 
with B'^ — E'^ — B'^/^'^. This condition is equivalent to the condition R^ Rf, where -R/ is 
the radius of the fast magnetosonic point (Beskin et al. 1998). 
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3. The flow near the axis 

Close enough to the axis, a jet may be formed if flow lines turn through even a small 
angle A^^ ~ 9. Let us consider the flow near the axis. To specify the conserving fluxes one 
should know the flow structure in the near zone. Note that at distances exceeding the 
flow fills all the space and plasma density and other parameters vary at the angular scale 
^ ~ 1. One can naturally assume that close enough to the axis, ^ <^ 1, both plasma density 
and the initial Lorentz factor are close to the values at the axis. Therefore we can take as 
a good approximation ry(^) = const and 70 — const. The flow lines of the pulsar wind are 
practically straight inside the fast magnetosonic point (but outside the light cylinder) and 
the flow is nearly quasiradial in this region. An example of such a flow is provided by the 
well-known split monopole solution (Michel 1973) 

* = vl>,^(l-cos^)«^*e.^^ B^^'M^^ E=-V*. (10) 

z Cri C 

Here it! and 9 are spherical radius and polar angle, correspondingly, and "^eq is the total 
(equatorial) poloidal flux. This solution is valid within the fast magnetosonic point. Taking 
into account that ^ oc at <C 1 in any quasiradial flow, one can describe any such a 
flow near the axis by an equivalent split monopole solution. Note that even in flows that 
are eventually collimated (e.g., Eichler 1993) the outflows well within the coUimation scale 
reduce to (10), which represents the precollimated state. 

The energy integral n should in principle be determined by the conditions at the fast 
magnetosonic point. Careful considerations show that deviations from quasiradial propa- 
gation remain small at this point in high a flows (Begelman & Li 1994, Tomimatsu 1994, 
Beskin et al 1998). Though these deviations are crusial for determining the position of this 
point, one can safely neglect them in the angular distribution of the parameters. Substituting 
equation (10) into equation (9), one gets 

n^^ea sin^^ 
= 70 H ■ 

Expressing back 9 via ^, we find the approximate energy integral as 

/i(^) = 70 {1 + ae, - (^/^e,)') } ~ 70 (1 + 2ae,*/^e,) , (H) 

where is the ratio of the Poynting to the matter energy flux at the equator. Because we 
are only interested in the flow near the axis, the "equatorial" values should be considered not 
as true values at the equator of the flow but as values at the equator of the split monopole 
flow that coincides close to the axis with the flow of interest. 
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Prom here on we will use the small angle approximation ip = ^ /'^eq <C 1, r <S z etc. 
We will consider the flow as ultrarelativistic, 7 1, because a Poynting dominated flow 
beyond the fast magnetosonic point is necessarily ultrarelativistic. Note that Poynting flux 
decreases toward the axis, the initial ratio of the Poynting flux to the kinetic energy flux 
being 



Thus, close enough to the axis, ip < (2(Teg)~^, the flow is matter dominated. This part of 
the flow is collimated already at the distance i? ~ i?j (Bogovalov 2001). We are interested 
in collimation of a wider flux tube where the flow is Poynting dominated. 

Note that decreasing of the Poynting flux toward the axis may be understood from 
a rather simple consideration that demonstrates that this feature is quite general and not 
a speciflc feature of Michel's solution. It follows immediately from equation (8) that the 
electric fleld decreases toward the axis. The toroidal magnetic fleld B^p oc J jrdr also 
decreases unless the current density j is singular at the axis. Therefore the Poynting flux 
EB^/Att decreases toward the axis. Our expression (11) for may be considered simply 
as expansion of general energy flux in small ip. 

The shape of the flux surfaces is conveniently described in cylindrical coordinates by 
the function r{il!, z) instead of -(/^(r, z). Then, e.g.. 



It follows immediately from Eqs (6, 8) that a oc rB/^ when 7]{ip) — const. Transforming 
derivatives to independent variables r and -0, one can write Eq.(5) as 



(12) 




In the small angle approximation the curvature radius may be presented as 




2 dr d'^r 



ahi(l + a) 



(13) 



9-0 dz'^ 



The energy equation (9) is now written, with the aid of Eq.(ll), as 



7 + 



riVV'l =7o(l + 2aegV')- 



V 



Taking into account that 7^ ^ 701^0 beyond the fast magnetosonic point (see, e.g., Beskin et 
al. 1998), one can substitute v = 1 into the denominator of the second term in the left-hand 
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side of this equation and write finally 



7 = 7o 



l + .,|2^-r(^) ' 



The system of equations is completed by the expression for a (sec Eqs.(6, 7, 12)): 



a — 



471777 



7o / dr 



7 



7 ydi/j 



(14) 



(15) 



4. Collimation of the flow 

Bogovalov (2001) very recently showed numerically that in the case of 7q < (Xeq the 
flow within the tube ip < {2creq)~^ is collimated into a cylindrical jet at the distance of the 
order of the radius of the fast magnetosonic point. One can see that within this flux tube 
(Jo < 1 and the fraction of the energy carried by the jet is only ~ a^^. This is insufficient to 
explain the observed jet in the Crab nebula. Here we consider whether a larger fraction of 
the poloidal flux (and, correspondingly, of the total energy) may coUimate at larger radii. 

The initial radial flow (10) may be presented near the axis as r = y/2ipz. Let us look 
for small deviations from this radial flow, 

r ^ ^/2^iJz{l - y); y('0,2;)-Cl. 

The approximate coUimation condition will be defined as follows. We shall solve the equa- 
tions retaining only terms of the lowest order in y and if we find y — > 1 for some poloidal 
field line, we conclude that this line coUimates. 



4.1. The case ^ < {2aeq)~^ 

It follows from Eq.(ll) that under the condition < (2(7eg)~^ the fiow is matter dom- 
inated. One can approximate this case by assuming cr <S 1, which is formally correct for 
ip <^ {2acq)^^. In this case the fiow evidently does not accelerate, 7 ~ 70, and a varies 
insignificantly, a ^ (Tq, as long as deviation of the flow from the purely radial one is small, 
y <C 1. Linearizing Eq.(13), one gets with the aid of Eq.(12) 



The solution to this equation is 



^In^. 

lo Rf 
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The initial point was chosen to he Rf because the equations are vahd only ai z ^ Rf and 
there is no coUimation a,t z <^ Rf (Begelman & Li 1994). One can see that if 

2^e. < JV, (16) 

where R is the radius of the termination shock, the flow is not coUimated at all. If the reverse 
is true, i.e. if 

2^e. > r^V, (17) 

the flow collimates at ip < {2aeq)^^. At 2o"eg > 7o the flow tube ip < {2aeq)^^ collimatcs 
already at 2; ~ Rf, which agrees with the results by Bogovalov (2001). At the condition (17) 
more flow than this may coUimate but as will be seen below, such additional coUimation 
invariably introduces acceleration, thus invalidating the approximation of the constant 7. 



4.2. The case t/j > {2aeq) ^, no acceleration 

In this region wc can take a ^ 1. Let us first assume that the fiow is not accelerated 
significantly and show that this implies no further coUimation. Substituting 7 = 70 and 
a — ao into Eq.(13) and linearizing it, one gets 



d'\zy)_ 
dz'^ 7o 



= 4. (18) 



The solution to this equation is 



In^ 



Because in the region of interest (Tq ^ 1 even a small deviation from the radial flow may 
result in significant growth of 7. The linearized energy equation (14) reads 



7 = 7o^l-4a.,V^^^j. (20) 

Our approximation 7 pa 70 is valid as long as the second term in the brackets remains small 
as compared with unity. Substituting the solution (19), one can see that our assumption of 
constant 7 is already violated at i? ~ i?/ if Aaeq > 7o; this case will be considered in the 
next subsection. By Eq.(20) the solution (19) is vahd at 

4^1nA<i. (21) 
7o 
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If it is valid until the termination shock, y remains less than unity for any ip > {2aeq)~^, so 
we come back to condition (16). If the reverse is true, the acceleration becomes significant 
at ^ 

zo^RfBx^pl^-^^ <R. (22) 

At this point y <^ 1. However, beyond this point coUimation may ensue. So let us consider 
coUimation in the accelerating flow. 



4.3. The case t/j > {2aeq) , accelerating flow 

Now one can neglect the first term in Eq.(20) and, substituting the Lorentz factor so 
obtained into Eqs.(13, 15) and neglecting the higher order terms in y, one obtains 

The solution to this equation with logarithmic accuracy (i.e. neglecting terms of [ln(i?/it!j)]~^) 
is 

((3/2) In 

y- 2V'(7o<7e,)V3 ■ (24^ 

Here Rq = Rf if 4o"eq > 7o and Rq = zq if (in-^j < Ao'eg/'yQ < 1. The fraction of the 
coUimated flux may be estimated from the condition y ~ 1 as 

((3/2) tof) 



4.4. Energy considerations 

Substitution of Eq.(24) into Eq.(20) yields 

/3 zV^^ 
7 = 2(^-7oae,ln— J , (26) 

which is consistent with the results by Tomimatsu (1994) and Beskin et al. (1998). In 
principle, the flow may accelerate up to the Lorentz factor 7mai;('0) — lo^o — '^loC'eqi^- At 
the coUimation point deflned from Eq.(24) as a point where y ~ 1, one obtains 7 ~ 7maa;('0)- 
So when a flux line turns through an angle ~ ^, a signiflcant fraction of the Poynting 
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flux is transferred to the plasma and there should be cr ~ 1 in the coUimated part of the 
flow. The collimation (and acceleration) process has enough time to be nearly completed 
only within the small fraction (25) of the poloidal flux. Outside this fiux tube the flow slowly 
accelerates according to Eq.(26). However, when the flow enters the termination shock, most 
of the energy still remains in the form of Poynting flux. 

Because Poynting flux decreases toward the axis of the flow, the fraction of the energy 
carried by the jet is small. Integrating the energy flux (11), one gets £'jet/-E'tot — ^-^ipc- 
the jet carries the fraction of the total energy 

^ 3 ((3/2) in A)-^' 



5. Discussion 

The fraction of the collimatcd flux is determined by the parameter 'jmax = ToCeg, which 
is Michel's (1969) magnetization parameter. This parameter is determined by the amount 
of plasma emitted by the pulsar. The density of the plasma generated within the pulsar 
magnetosphere is conveniently normalized by the Goldreich- Julian charge density: 

Kpoj kBD, 
e 27rec ' 

where B is the magnetic field in the magnetosphere, k the multiplicity factor. Within the 
light cylinder the field is nearly dipolar, B — /i/R^, where /i is the magnetic moment of the 
star, and conservation of the mass flux imphes n oc B. Beyond the hght cylinder n oc 1/B?. 
The Poynting flux in the pulsar wind may be expressed via the magnetic field at the light 
cylinder, Bl — jiVt^ /c^, as 

Now one can estimate the magnetization parameter as 

^^^^ = ^ = ^ 3 X 10^-^ (28) 

where = lO^^^^o G-cm^, k = 10'^/«4 and P is in seconds. For the Crab pulsar, where 
P — 0.033 s and /^ao = 5, we have 

_ 6 • 10^ 

^max ■ V / 

Prom the observed amount of radio emitting electrons (those with the energy > lOOMeV) 
one can get an upper limit for the injection rate of these electrons about ~ 10^^ s~^, which 
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corresponds to k ~ 3 • 10^. If the pulsar ejects a larger amount of pairs, their energy in the 
nebula should be less than 100 MeV. 

The standard models of plasma production in pulsars (see, e.g., the review by Arons 
1983) typically predict k ~ 10^ — 10^ and 70 ~ 100. Different modifications were proposed 
to explain the observed powerful gamma radiation (Cheng, Ho & Ruderman 1986, Usov & 
Melrose 1996, Lyubarskii 1996); these models predict significantly larger plasma production. 

However, we are not aware of models predicting k, > few x 10^. If protons are also emitted 
by the pulsar, the corresponding magnetization parameter is reduced by a factor m/nip. 
However theoretical arguments (Arons 1983) favor ~ 1. 

It is very difficult to place limits on the wind parameters directly from observations. 
Limitations found by Wilson & Rees (1978) who studied induced scattering of the radio 
emission in the pulsar wind were obtained on the assumption that the cyclotron frequency 
near the light cylinder is less than the radiation frequency. However the reverse is true in the 
case of the Crab pulsar (Barnard 1986, Lyubarskii & Petrova 1996). The scattering rate in 
the magnetized plasma is small if the pulsar radiates only waves polarized perpendicularly 
to the local magnetic field (Blandford & Scharlemann 1976, Lyubarskii & Petrova 1996). In 
this case the scattering begins only beyond the polarization limiting radius, which lies, in 
the Crab case, well outside the fight cylinder. 

Because pairs are produced by gamma radiation, one can reasonably conclude that the 
pair luminosity does not exceed the gamma luminosity. This places the limit Ueq > 100. One 
can see from Eq.(27) that Ejet/Etot ~ 10~^ even at Oeq — 100 and 70 = 1. 

Although the fraction of total luminosity that emanates from the jet is not yet published, 
and in any case is not quite the same as the fraction of mechanical output that emerges in the 
jet, visual inspection seems to imply that 10~^ is a conservative guess. The implication that 
(Teq is at most 10^ and 70 ~ 1 runs counter to conventional pulsar wisdom. Probably observed 
constraints on relativistic beaming will provide an even sharper constraint on and 70, 
making the observed jet even more confusing for conventional pulsar theory. However, since 
existing theory for pulsar magnetospheres is not rigorous, perhaps it needs to be reexamined. 
The jet may be formed via hoop stress collimation of the the pulsar wind only if the wind 
is much more mass-loaded than previously believed. The alternative possibifity, that the jet 
is formed already near or within the pulsar magnetosphere, R Rl, also requires serious 
modifications to the conventional picture of the pulsar magnetosphere. 

We thank the Israeli Science Foundation and the Arnow Chair of Physics for their 
generous support. 
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